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ABSTRACT 

While second and higher order correlations of the light distribution have re- 
ceived extensive study, the lowest order probability distribution function (PDF) 

- the probability that a unit volume of space will emit a given amount of light 

- has received very little attention. We estimate this function with the aid of 
hydrodynamic simulations of the ACDM model, finding it significantly different 
from the mass density PDF, and not simply related to it by linear bias or any 
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of the other prescriptions commonly adopted. If the optical light PDF is, in 
reality, similar to what we find in the simulations, then some measures of Q m 
based on mass-to-light ratio and the cosmic virial theorem will have significantly 
underestimated Q m . Basically, the problem is one of selection bias, with galaxy 
forming regions being unrepresentative of the dark matter distribution in a way 
not described by linear bias. Knowledge of the optical PDF and the plausible 
assumption of a log-normal distribution for the matter PDF will allow one to 
correct for these selection biases. We find that this correction (which amounts 
to 20 — 30%) brings the values of Vl m estimated by using the mass-to-light ratio 
and the cosmic virial theorem to the range fl m = 0.2 — 0.3, in better agreement 
with the WMAP result than the uncorrected estimates. In addition, the rela- 
tion between mass and light PDFs gives us insight concerning the nature of the 
void phenomenon. In particular our simulation indicates that 20% of mass is 
distributed in voids, which occupy 85% of volume in the universe. 

Subject headings: cosmology: large-scale structure of Universe - - cosmology: 
theory — galaxies: formation — methods: numerical 



1. Introduction 

In modern cosmological studies, the distribution of the matter density has received 
extensive study, based on theory and on increasingly accurate dark matter numerical simu- 
lations. We may study the matter distributions directly or by looking at the distribution of 
halos after applying group finding algorithms to identify virialized "halos" in N body simu- 
lations. The simplest statistic is the probability of a sphere of radius R to contain mass in 
the range M — (M + dM) or, alternatively, to have a density in the range p m — {p m + dp m ), 
dividing by the volume of the sphere. We may more conveniently use the normalized dimen- 
sionless density, y = p m / (p m ), as the parameter, further dividing by the mean mass density, 
and consider the probability distribution function (PDF) of y. By the definition of a PDF, 

POO 

/ f R (y)dy=l, (1) 
Jo 

and the commonly used density variance is 

POO 

ol= (y- l) 2 f R (y)dy = ((y - lf) R = (S 2 m ) R , (2) 
Jo 

where the additional definitions are given to establish nomenclature used in this paper, and 
the subscript R shows the dependence of the averaged quantities on scale R (we usually drop 
R, however, unless we need to emphasize it). 
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It is known in CDM models that the logarithm of the density shows a normal distribution 
to a good approximation when the structure formation has reached the non-linear regime 
(Coles & Jones 1991; Kofman et al. 1994; Taylor & Watts 2000; Kayo, Taruya, & Suto 2001): 



The parameter that specifies the distribution, u, is related to the variance a of the density 
field as, 



It is known that this log-normal function can be obtained from the one-to-one mapping 
between the linear random Gaussian field and the nonlinear density field (Coles & Jones 
1991), although the physical meaning of the transformation is not well understood. The 
mapping, nevertheless, seems to capture an important aspect of the nonlinear evolution of 
the density field in the universe. 

In parallel to the PDF of the mass distribution f(y), we also define the PDF of the light 
distribution by g(j), where j is the normalized luminosity density. It is often convenient to 
define the logarithm of each parameter as Y = logy and X = log j. 

It is somewhat curious that much less attention has been paid to the PDF of the observed 
light in the universe, whereas two particle and higher order distribution functions have been 
studied extensively in the literature. The work on "counts-in-cells" by Efstathiou et al. 
(1990) implicitly deals with the projected mass PDF in the linear regime, and the recent 
work of Dekel & Lahav (1999) studied this function in the context of "stochastic biasing". 
The void distribution function, which depends on all higher order distributions, has received 
some attention (e.g. Vogeley, Geller, & Huchra 1991; Vogeley, et al. 1994; Benson, et al. 
2003). However, we are not aware of any systematic analysis of the optical PDF or attempts 
to relate it with the mass PDF (but see the work of Sigad, Branchini, & Dekel (2000) 
where they attempted to relate the galaxy biasing function and the cumulative distribution 
functions). 

In this paper, we consider the relation between g(j) and f(y), or more explicitly, the 
relation between X and Y, X — X{Y). These relations can be determined by the physics 
of galaxy formation, and we illustrate them with the use of hydrodynamic simulations of 
the concordance ACDM model. The results allow us to elucidate the nature of the void 
phenomenon (Kirshner et al. 1981; Rood 1988 and references therein; Peebles 2001) in the 
CDM model that is confronted with the observation. Our results indicate that the relation 
between X and Y is not linear, contrary to what has usually been assumed. This also points 
to the presence of systematic biases, when galaxies are used as tracers of the mass, in the 




(3) 



c4 = ln(l + 4). 



(4) 
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conventional estimates for the total mass density of the universe and for the evaluation of 
the cosmic virial theorem. Our results may be used to correct for these biases. 

We show in Section 2 the PDF of the mass and light distributions and their relation 
to each other using a hydrodynamic simulation. In Section 3 we discuss implications to 
observations, in particular, how selection bias can be corrected with the knowledge of these 
functions. 

2. PDF from Hydrodynamic Simulations 
2.1. Simulation 

The simulation we use is similar to Cen & Ostriker (1993), but has been significantly 
improved over the years by adding more elaborate modeling. The simulation which uses the 
Total Variation Diminishing (TVD) method (Ryu et al. 1993) is performed for a box-size of 
Lbox = 25/i _1 Mpc with Eulerian hydrodynamic mesh of 768 3 and cosmological parameters 
of (fi m , tt A , tt b , h, a s ) = (0.3, 0.7, 0.035, 0.67, 0.9), which are close to the estimate by WMAP 
(Spergel et al. 2003). This is the same simulation that has been extensively used in our 
earlier papers (Nagamine, Fukugita, Cen, & Ostriker 2001a,b; Nagamine 2002) to discuss 
star formation history, luminosity function, and the nature of Lyman-break galaxies. We 
refer to Cen et al. (2002) and Nagamine et al. (2001b) for more details of the simulation. As 
in our previous analyses we employ the population synthesis model of GISSEL99 (Bruzual 
& Chariot 1993; Chariot 1999, private communication) to calculate the stellar luminosity. 
Unless otherwise noted, all the luminosity used in this paper is in the B band, neglecting 
the dust extinction effect. 



2.2. One-Point PDF of Mass and Light 

We calculate mass and light contained in spheres of radius R located at arbitrary loca- 
tions, i.e., smoothed with a tophat filter. We then compute the probability dP of a sphere 
to have the normalized density parameter y — (y + dy) for mass and j — (j + dj) for 
light: dP = f(y)dy = g(j)dj, where f(y) and g(j) are the one-point PDF of mass and light, 
respectively. 

In Figure 1, we show the one-point PDFs f(y) in panel (a) and g(j) in panel (b) 
computed from our simulation with R = \h~ x Mpc, a scale which is small enough to be 
sensitive to the local galaxy formation environment but large enough for the simulation to 
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be accurately representing physical values. Since our binning is done in logy, the plotted 
histogram actually corresponds to (In 10) yf(y) [f(y)dy = (In 10) yf(y)d\ogy], and similarly 
for g(j). Note the difference in the scales of the abscissas of the two panels. 

The solid curves plotted over the histogram in both panels are log-normal distribution 
with the variance parameter u 2 computed from the same simulation, u 2 = 3.1 for dark 
matter and u 2 = 5.6 for light. The agreement between the computed dark matter PDF and 
the log-normal distribution is reasonably good. The departure in our computation, especially 
at small logy, is ascribed to the limited box-size of our simulation. Note that the light PDF 
strongly deviates from the log-normal distribution. (The long- dashed line in Figure 1(b) will 
be explained in Section 2.3.) 



2.3. Relation between logy and logj 

Figure 2 shows the light density parameter X = logj against the dark matter density 
parameter Y = logy. 1 The point in the figure represents mass and light contained in each 
sphere. The solid line shows the median value in each X bin, with error bars in X-direction 
being the quartiles on both sides. The long-short dashed line shows the identical regression 
X — Y, and the long-dashed line shows an empirical fit to the simulation data: 

X = aY ~ JV^cj~ 2 ' {Y>C) - (5) 

The parameter c corresponds to an effective cutoff of the mass density, below which the 
sphere does not contain any galaxies that emit light. The figure indicates that c ~ —0.5, i.e., 
the region of mass density smaller than 1/3 of the mean is dark. This dark region occupies 
approximately 60% of the volume when smoothed with a tophat sphere of R = \h~ x Mpc, 
and contains a fraction ~ 8% of the total mass (these fractions are larger for voids that are 
operationally defined in observations: these numbers will be derived in Section 2.4). This is 
identified as the "void phenomenon" in the CDM universe. Luminous galaxies reside only 
in the dense regions. Sub-luminous galaxies prefer to lie in relatively low density regions, 
but the sharp increase of logj at around Y m c means that they live only in the edge of 
the voids, certainly not in the middle of the voids. These trends seem to be qualitatively 
consistent with the observed void phenomenon (Peebles 2001). 



: We note that a similar figure was reported by Blanton et al. (1999) from our earlier simulations, but 
they used stellar mass density instead of j. The figure was also truncated at (logy, logj) w (0.5,0.5). The 
simulation box-size we use here is smaller than that Blanton et al. used, and does not contain massive 
clusters that show the suppression of luminosity density at the highest mass density region. 
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Given the relation between X and Y in Equation (5), and assuming a log-normal dis- 
tribution for f(y) which is specified by ou, one can compute the PDF of light via g(j) = 
f(y)dy/dj. The long-dashed curve given in Figure 2 is a fit by a = 1.4, 6 = 0.7, and 
c = —0.85 for the range —0.85 < Y < 2.5 (the actual value of c is smaller than —0.5 because 
of the slow asymptotic nature of the curve near the critical value). The result of g(j) thus 
obtained is shown as the long- dashed line in Figure 1 (b). The agreement between this model 
calculation and the simulation is good, verifying the self-consistency between the two PDFs 
in Figure 1 and the relation between {X,Y) given in Equation (5). 

We note that there are two constraints on g(j), 

POO POO 

/ 9ti)dj = l, / .///(./>// i. (6) 

Jo Jo 

The empirical fit we adopted in Equation (5) alone does not satisfy the first constraint of 
(6); we must add the zero luminosity component, i.e., a Dirac delta function S D (j). Hence 
the function g(j) should be written as 

9(j)=A5 D (j)+g'(j), (7) 

and the first constraint reads 

poo 

A + / g'(j)dj = 1, (8) 
Jo 

where g'(j) corresponds to the light-emitting component discussed above, i.e., the long- 
dashed line or the histogram itself in Figure 1 (b). We find A = 0.6 (see Section 2.4) and 
that the two constraints are satisfied for the model curve within an accuracy of ~ 5% where 
we limit the integral to log j < 2.5, the shot-noise limit of our simulation. Needless to say, 
the simulation data automatically satisfy the constraints by the definition of PDFs. 



2.4. Distribution of Matter and Light 

In this subsection, we study the distribution of mass and light from different angles and 
consider the relation to observations. Figure 3(a,b) shows the cumulative volume fraction 
of regions with mass density < logy and with light density < logj, respectively; i.e., the 
probability [P(< logy, logj)] that a volume element has mass (light) density less than a 
given value of logy (logj) in the abscissa. The solid line is computed from the simulated 
PDF, and the short-dashed line is obtained by integrating the model curve of Figure 1. We 
have seen that no luminous objects are present in the region of logy < —0.5. 
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Figure 3(a) shows that few regions are emptier than logy < —2: i.e., the "voids", while 
they are completely dark, are not thoroughly devoided in mass (cf. Peebles 2001). From this 
figure we see that the region with logy < —0.5 occupies 60% of the volume (dotted line). 

We see in Figure 3(b) that the light density exhibits a smooth distribution until it 
shows a jump at the faint end. The magnitude of the jump at log j = —6 (solid line) means 
that spheres in 60% of volume contains no light; hence the value of A in Equation (7) is 
determined to be A = 0.6. The lower limit of logj is an artifact of the resolution of the 
simulation, but we consider that the realistic cutoff as determined by the Jeans and cooling 
conditions is rather close to this limit. The flattened distribution of P(< logj) towards the 
low luminosity end indicates that such a jump is present in reality. 

Let us consider voids defined in observations. Most surveys sample galaxies down to 
about 4 mag below L*. By integrating a standard Schechter luminosity function (approxi- 
mately 0* = 0.02/i 3 Mpc -3 , M* = -20 + 51og/i, and a = -1.25) from 0.025L* (i.e. magni- 
tude limit of M — —16) to infinity, we expect that a sphere of radius R = 1 Mpc contains 
0.4 galaxies (here we assume h — 1 for simplicity). This means that a sphere of luminosity 
< 0.01L* cannot be observed in current standard galaxy surveys. Since the luminosity in a 
sphere of R = 1 Mpc with the above Schechter function (with a magnitude limit of 0.025L*) 
is 0.08L*, the observational threshold of light density parameter is j t h = (0.01/0.08) = 0.125, 
or logyth ~ using Figure 2. Namely, spheres in the region of logy < are observed as 
'dark'. Using Figure 3(a) this means that 85% of volume is identified as voids (dotted line). 

In order to see the mass in voids, we plot in Figure 4 (a) the mass fraction and (b) light 
fraction that are contained in regions of < logy, i.e., f£ y'f(y')dy' and Jq j'g(j')df for dark 
matter and light, respectively. Again, the solid line is computed from the simulated PDF, 
and the short-dashed line is from the model curve. The result is robust against the change 
in the lower limit of the integration. The discrepancy between the short-dashed line and the 
solid line is due to the difference in the shape of the simulated and the model PDFs. 

From Figure 4(a) we read that 8% of mass (corresponding to logy < —0.5 where the 
sharp cutoff was observed in Figure 1 ) does not emit light at all, and observationally defined 
voids (logy < 0) contain 20% of the total mass as indicated by the dotted lines. The relation 
between the mass fraction and the light fraction is represented more clearly in Figure 4(c), 
which indicates that 90% of light comes from higher density regions where 50% of mass is 
contained (or in other words, only 10% of light is contained in lower density region where 
50% of mass is contained; see dotted lines). The departure from the diagonal line in panel 
(c) shows the difference in the manner that dark matter and light are distributed in the 
universe. 



- 8- 



3. Applications 

3.1. Mass-to-Light Ratio 

The result we have seen in the preceding section indicates that the conventional evalu- 
ation of the the mass density from the luminosity density (p L ) and the mean mass-to-light 
ratio (M/L), as (p rn ) = (pl) (M/L) (Peebles 1971), underestimates the true value. 

If galaxy formation is considered as largely a local process, then it is expected to be 
primarily a function of local gas density and gas temperature. To a first approximation 
baryonic gas density should roughly trace the total mass density, and indeed this is verified 
in our hydrodynamic simulation as shown in Figure 5. Therefore, let us assume that galaxy 
formation is primarily a function of only the total mass density p. Namely, we take the 
normalized light density j and the local mass-to-light ratio p{y) = ^(y) = py/j(y) as 
functions of y. 

The volume averaged mean mass density can be calculated by 

P oo 

(p) v = / pyf(y)dy, (9) 

Jo 

and the volume averaged mean light density is 



/■oo 

(j) v = / j(y)f(y)dy- (10) 

JO 

So the appropriate global mean mass-to-light ratio would be 

f QO pyf(y)dy 

Wv ~ J °°j(y)f(yW { } 

which is defined so that, by construction, 

(P)v = (P)v U)v ■ (12) 

In observations one cannot observe the dark matter directly, and therefore, one samples 
only where the light is. This could lead to an underestimate of mass since volume elements 
in the universe are not treated equally. If we evaluate (p) L observationally, (p) Lohs , by 
studying only those regions which emit light and weight them by the emitted light j, then 
we are effectively evaluating the integral 

, * J °°pyjf(y)dy 

Jo 00 pyjgiM (14) 

f °°f9U)dj ' 
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To see how much underestimate could be caused by weighting by light, take the following 
simple model in which we assume that the mass-to-light ratio is constant (= /xq) in the 
galaxy forming regions. Namely we assume that the light density is linearly proportional to 
the mass density, but that there are empty voids: 

j — for y < y 1 ( "voids" ) 

j = —y for y > y\ ( "galaxy forming regions" ) . (15) 

This toy model is shown as short-dashed line in Figure 2. The mean mass-to-light ratio we 
defined in Eq.(ll) becomes 

{ " )v = ""l^JWv (16) 
(, , Io 1 yf(y)dy \ , 

= »[ 1+ ttjm)- (17) 

The luminosity weighted average [Eq.(13)] then becomes simply 

So the true volume weighted average {fj) v is larger than the luminosity weighted average 

(^)L,obs b y a factor F: 

r . My 1 J^yf(y)dy 

The ratio of the integrals in Eq.(19) is just the mass fraction in "voids" divided by that in 
"galaxy forming regions". We estimated that the "voids" contain 20% of the total mass; this 
means that the correction factor will be T = 1.25. This leads to an underestimate in Q m by 
the same factor. 

The estimate of Vt m obtained by using (p B ) = (2.2±0.4) x 10 8 hL Q>B Mpc~ 3 and (M/L) = 
290 ± 35/i (McKay et al. 2001) from gravitational lensing shears, Vl m = 0.23 ± 0.05, should 
therefore be corrected to give Vl m = 0.29 ± 0.06, in good agreement with the value from the 
WMAP experiment: Vl m = 0.29 ±0.07 (Spergel et al. 2003). The same correction also applies 
to similar estimates done to date, which always show rather small values of Q m compared to 
the estimate of WMAP. For example an estimate of Bahcall et al. (2000), Q m = 0.16 ± 0.05, 
using mass-to-light ratio as a function of scale is modified to 0.20 ±0.06 after the correction. 
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Similarly a low value of Vt m = 0.15 obtained by Fukugita, Hogan, & Peebles (1998, equation 
(28)) based on galaxy luminosities should also receive an upward 25% correction. 

If the light PDF g(j) is estimated directly from the observation, one can numerically 
integrate the relation dj/dy = f(y)/g(j) to obtain j(y) using a cosmological simulation. With 
j(y) one can easily obtain the mass-to-light function /i(y) = py/j(y) and the correction factor 

r (vh ( f™yf(y) d y \ , ( J™yif(y) d y \ , 9m 
Ml \Io 00 j(y)f(y)dy) 1 K^FfWy) ' [ } 



3.2. Cosmic Virial Theorem 

Let us now turn to the application to the "cosmic virial theorem". Peebles (1980, 1993) 
shows that, if one assumes that galaxies trace mass, the "cosmic virial theorem" which relates 
the relative velocity dispersion and the mean mass density can be approximated by 

2(3- 7 ) a{rf 

n - ~ 3 Wy^- (21) 

To derive this equation, (, gg (r) = i gm (j') is assumed, where (, gg (r) = (r/r 0tgg )~' y is the galaxy- 
galaxy correlation function, and £ gm (r) = (r/r 0tgm )~ 13 is the galaxy-mass cross-correlation 
function. 

As we showed earlier, galaxies do not trace mass well. Therefore if instead one were to 
use the galaxy-mass cross-correlation function in deriving the above relation, the value of 
Q m would have been larger or smaller by a factor of 

-,/ ^rn,grn A 3 @\ ( ^0,gg 



TV _ " "m,gm _ i ~ r_ u, gg /<-,<-, x 

where r is measured in units of hr l Mpc and the right-hand-side of Eq. (21) is measured at 
r = l/i _1 Mpc. 

The values of 7 and r 0>gg are observationally known (7 ~ 1.8, r , gg ^ 5/i _1 Mpc) and (3 
and r 0i9m can be estimated in the linear regime as follows: 
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where we have assumed a relation 



and therefore, 



(pg) \ (Pm) 



(25) 
(26) 

(27) 



5A 4 5 A (28) 



P ) 9 \P 

The analytic fit that we adopted in § 2.3 Eq. (5) implies k ~ a ~ 1.4 for the domain of 
< logy < 2 that concerns us. By comparing Eq. (26) and the original functional form for 
£ 9m , one obtains 

r , gm = k~ r 0>gg . (29) 
For any value of k larger than unity will therefore gives 

P = 7 and r , 9m < r , Sfl . (30) 

Using these relations in Eq. (22), the matter density of the universe estimated by Peebles 
by this method (Q m ~ 0.2) should be increased by a factor of V = 1.2, which yields a 
corrected value of Q m = 0.24. The agreement with the WMAP result is improved by this 
correction. 



4. Discussion and Conclusions 

Very little work has been done on the optical probability distribution function (PDF): 
the probability of finding a region (of some given size) containing a specified range of optical 
light output from galaxies within that region. Even less has been done in relating the optical 
PDF with the matter (or dark matter) PDF. 

What is known is that the low density regions (optical "voids") are relatively unpopu- 
lated by galaxies (Kirshner et al. 1981; Rood 1988), and consequently are optically faint as 
compared from the dark matter simulations of the currently popular concordance model. If 
this is correct, then "bias" — the ratio of galaxies to dark matter — is different in high and 
low density regions. And more seriously, our observational tracers of cosmic structure are 
giving us a view of the universe that is inherently contaminated by selection bias. 
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We compute the effects in a TVD cosmic hydrodynamic simulation and then estimate the 
error that would have been involved, had these simulations been observed and analyzed using 
conventional techniques. Focusing on the global value of Q m (which of course is known for the 
simulated universe), we can estimate the corrections, and when applying these corrections 
we obtain increased values of fl m to the range 0.2 — 0.3, which is in better accord with the 
CMB determinations. 

The relation between mass PDF and light PDF also gives us insight concerning the 
nature of the void phenomenon, and enables us to estimate the mass and volume fractions 
occupied by voids. Our result indicates that no galaxies can be present in the middle of voids: 
sub-luminous galaxies are distributed in the edge of voids, whereas luminous galaxies reside 
only in high-density filaments. The quantitative basis for these statements will be presented 
in a subsequent paper, but the concept is as follows. A large void can be considered as a 
small piece of a lower density universe. The implies a shift to a lower characteristic halo 
mass in such regions and a shift in the horizontal scale in Figures (2) and (3a) to domains in 
which galaxy formation is less efficient. At the centers of such voids all galaxy formation is 
suppressed and at the edges high luminosity galaxies (from higher mass halos) are suppressed 
more than lower luminosity systems. A recent numerical study by Gottloeber et al. (2003) 
supports this view. 

Future work using large-scale surveys can provide the empirical data needed to compute 
the optical PDF properly. When this is combined with detailed numerical simulations of the 
dark matter, we will be able to self-consistently evaluate these effects better and use the so 
determined relation between light density and mass density in our analysis of observational 
data. 

We thank Michael Strauss for useful comments on the manuscript. K.N. thanks for the 
hospitality of the Institute of Astronomy, University of Cambridge where a part of this work 
was completed. This work was supported in part by grants AST 98-03137 and ASC 97-40300 
in Princeton, and Grant in Aid of the Ministry of Education 13640265 in Japan. 
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Fig. 1. — One-point PDFs of dark matter f(y) (panel (a)) and of light g(j) (panel (b)) 
computed from a cosmological simulation of box-size 25/i _1 Mpc with a tophat window of 
size R = l/i^Mpc. Note that the plotted histogram actually corresponds to (In 10)yf(y) 
because the binning is done in log(y) instead of linear y. The solid curve in both panels 
is the log-normal fit with the variance parameter uo 2 = 3.1 (for dark matter) and 5.6 (for 
light) [see eqs. (3) and (4)] computed from the simulation. The long-dashed line in the lower 
panel is computed by assuming the log-normal distribution for f(y) and the analytic fit to 
the relation between logy and logj as shown in Figure 2 and described in § 2.3. 
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Y=log(y) 

Fig. 2. — Relation between X = logj and Y = logy, where j and y are the normalized light 
and matter density parameter, respectively. Each point in the figure represents a tophat 
sphere of radius R = lh^ 1 Mpc placed in the simulation box arbitrarily. Notice that in the 
significant fraction of the volume having mass density less than 1/3 of the average value, 
there is effectively no light density, i.e., M/L — ► oo. Solid line: median values in each X 
bin, and the error bars are the quartiles on each side. Long-short dashed line: linear relation 
X — Y. Long dashed line: empirical fit explained in § 2.3. Short-dashed line: toy model 
described in § 3.1 with a cutoff at logy = —0.5. 
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Fig. 3. — Volume fraction of regions with mass density < logy (panel (a)) and with light 
density < logj (panel (b)); i.e, the probability that a volume element has mass (light) 
density less than a given value of logy (logj) in the abscissa. The solid line is computed 
from the simulated PDF, and the short-dashed line is obtained by integrating the model 
fitted PDF shown in Figure 1. The two dotted lines in panel (a) indicate that the region 
with logy < —0.5 (logy < 0) occupies 60% (85%) of the volume. 
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Fig. 4. — Mass fraction (panel (a)) and light fraction (panel (b)) that is contained in regions 
of < logy. Note that the light is more concentrated in high density regions than the dark 
matter is. The solid line is computed from the simulated PDF, and the short-dashed line 
is obtained by integrating the model fitted PDFs: J y'f(y')dy' and f£ j'g(j')df for dark 
matter and light, respectively. The two dotted lines in panel (a) indicate that the region 
with logy < —0.5 (logy < 0) contain 10% (20%) of the total mass. Panel (c) shows the mass 
and light fraction for the same values of logy as given in panels (a) & (b). The deviation 
from the diagonal long- dashed line indicates the difference in the manner that dark matter 
and light are distributed. The two dotted lines in panel (c) indicate that 90% of light comes 
from higher density regions where 50% of mass is contained (or in other words, only 10% of 
light is contained in lower density region where 50% of mass is contained). 
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Y=log(y) 

Fig. 5. — Relation between gas overdensity p gas / (p gas ) and total mass overdensity Y = logy, 
where y is the normalized matter density parameter. Each point in the figure represents a 
spherical volume element of radius lh" 1 Mpc in the simulation box. The solid line is the 
median values in each Y bin, and the error bars are the quartiles on each side. The diagonal 
long dashed line shows the linear relation between the two quantities plotted. 



